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DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD 
FOR PARABOLIC EQUATIONS 

H. Kaneko, K. S. Bey and G. J. W. Hou 
Abstract 

In this paper, we develop a time and its corresponding spatial discretization scheme, 
based upon the assumption of a certain weak singularity of ||wt(t)||^ 2 ( n ) = ||u f H 2 , for the dis- 
continuous Galerkin finite element method for one- dimensional parabolic problems. Optimal 
convergence rates in both time and spatial variables are obtained. A discussion of automatic 
time-step control method is also included. 
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1 Introduction 

In this paper, the following standard model problem of parabolic type is considered: 

Find u such that 

u t (x,t) — u xx (x,t) = /(x, £), x G Q, t G R + , 

u(x,t) ~ 0, x G G R + , (1*1) 

u(x : 0) — uo(x), x G ft, 

where ft is a closed and bounded set in R with boundary <9ft, = (0,oo), u xx ~ d 2 u/dx 2 , 

u t = du/dt, and the functions f and uq are given data. 

In this paper, we take advantage of a certain regularity of \\utW 2 , to discretize the time 
variable in a way that the optimal convergence rate of the discontinuous Galerkin (DG) method 
is obtained. In a series of papers [1], [2], [3] and [4], Eriksson and Johnson presented adaptive 
DG finite element methods in which time and space variables are adjusted using a posteriori 
estimates of the error. The current method is concerned with establishing a priori error estimate 
for the DG method which is optimal. 

The following conditions will be assumed. Let (/i, T, S) denote a finite element discretization 
satisfying 
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1. h is a positive function in C l (ft) such that 

I V h{ x )\ < M, for all x € fl and for some M > 0. 


2. T = {if} is a set of triangular subdomain of ft with each triangular element having 
diameter hx such that 

c\h\< [ dK for all K e T, 

Jk 

and associated with the function h through 


C 2 hx < h(x) < hx, for all x e K, K € T , 


where c\ > 0, C 2 > 0. 

3. S is the set of all continuous functions on ft which are polynomials of order r in x — (zi, X 2 ) 
on each K E T and vanish on dft. 

For the DG method for (1.1), we partition R + as 0 = to < t\ < • * * < t n < • • - where we let 
I n = (t n -i,t n ] with k n ~t n — t n ~ 1 . For each time interval, with q a nonnegative integer, define 

W hk = {v: R+ - V m : v\ In € P g (I n ), n = 1, . . . Nj, 


where 



the space of piecewise linear splines 
defined over ft = [0, L\ with breakpoints 
0 — xo < x\ < * - < x m = L 
and h m — max()< 1 • j 


and 

Q 

Pq(In) = Mt) = ^2 v it l:v i € V m }. 
i=0 

Even though linear splines are used in the definition of V m , splines of any order can be used to 
form V m . 

The DG method is defined as follows: 


Find U such that for n = 1 , 2 , • * • , U\axi n € Whk and 

/ {(U t ,v) + a(U,v)}dt + ([U} n ^ 1 ,v+_ 1 )= / (f,v)dt for all t, € (1.2) 

Jin Jin 
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where [w] n = iv„ - w~, ( * = lim s _ +0+( _ ) w(t n + 5 ), 1/(7 = u 0 , (u,v) = f Q u(x)v(x)dx and 

a(u, u) = (VC/, Vu). 

As stated earlier, the approach here is to first estimate ||it t || 2 . For example, if uq(x) = tt — x 
in (1.1) with f(x, t) = 0 and = (0, 7r), then the actual solution is given by 

00 

u(x, t)=^2 u°e~ j2t sin (jx), (1.3) 

3 = 1 

where 

m“ = f Jo " (tt - a:) sin(jx)dx 
= l{f-j 2 srn 2 jn} 

= 0 ( 1 ). 

In the following, C’s denote generic constants whose values change as they appear. From (1.3) 

OO 7 DO 

IMOIll = IKWIlLw = J2cfe- 2ft = J t E^" 2j2t - (1-4) 

j=l j=l 

The last equality in (1.4) is justified because of the uniform convergence in t of YljLi Ce~ 2j2t . 
Now using the fact that f£° e~ x2 dx < oo, a simple change of variables (say, y = jy/2t) will show 
that the last expression in (1.4) is which leads to 

lk(t)l| 2 = o(*- 3/4 ). 

A similar argument shows that if Uj = O(jz) for some initial value function uo(x), then 
\\u t (t)\\ 2 = 0(t -1 / 4 ). This case arises when uq(x) — min(x,7r — x) for x £ (0,7r). If u j de- 
cays faster than j~ 2 5 as j — » oo, then ||ut(i )||2 will be bounded as t 0, and in this case, time 
interval can be partitioned uniformly. Before we present the current method in the next section, 
one of the main theorems (theorem 2.3) from [2], which is pertinent to this paper, is recalled 
below. 

Theorem 1.1 (Eriksson and Johnson [2]) Let u be the solution of (LI) and U that of (1.2). 
Assume that V m C V m ~i and k n < 7 & n +i for all n and for some 7 > 0. Then there exists 
constants C only depending on c\ and C 2 (see condition 2 above) such that for q = 0,1, and 
N = 1,2,..., 

11“ - U ||/„ < CL n max E qn (u), 

l<n<N 
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and for q = 1, N — 1,2, . . . , 


11“ (t N ) - U N || 2 < CL n max . £ 2 „(u), 

l<n<N 

where Ljv = (Iog(t^/fe^) + l) 1 / 2 , 

£?n = min fc^||^t J) ||/ n + \\hlD 2 u\\ In , 9 = 0,1,2 (1.5) 

J<9+1 

with u[ l) = u*, = Au« and \\w\\j n = max* €/n ||u;(t)||2- 

The term minj< g+ i fc£||u^||/ n in (1.5) describes the error associated with discretization in 
time. Thus, if ||«^||/ n is bounded for each n and j = 1,2,3, then the DG method is of the 
jth order accuracy in time. In many cases where the initial condition is incompatible with the 
boundary conditions, \\u t \\j Y is unbounded as t — * 0 + . We demonstrate in Section 2 that the 
graded time partitions can be established to restore the optimal rates of convergence even in 
the presence of non-smooth data. It should be pointed out that similar non-uniform graded 
discretization schemes were used in capturing the solutions of a class of weakly singular integral 
equations, -e.g., [8], in capturing the solution of parabolic integro-differential equations with 
memory term, -e.g., [9], in capturing the solution of a class of initial value problems, -e.g., [13] 
as well as in capturing the solution of linear parabolic equations in the /i- vers ion of the DG 
finite element method, -e.g., [12]. Some other notable papers in this topic are, -e.g., [5], [10] and 
references cited therein. The recent paper by Schotzau and Schwab [12] is of particular interest 
here, since it also deals with the /i-version of the DG method for parabolic problems with 
non-smooth data, particularly the topic covered in section 5.2 of [12]. Schotzau and Schwab 
derive a graded time partition scheme by analyzing a function defining the initial condition. 
More specifically, uq in (1.1) is examined by the K- method of interpolation as a function which 
belongs to an intermediate space between i/ 2 , the Sobolev space, and L 2 and the graded time 
partitions are chosen accordingly. In this paper, graded time partitions are derived in Section 2 
by analyzing the regularity of ||u t || 2 . The difference between the current approach and that of 
Schotzau and Schwab in [12] manifests itself in an example (Section 3) in which we demonstrate 
that the current method gives rise to more relaxed time partitions. The paper [12] by Schotzau 
and Schwab is a must-read for anyone who is interested in the DG finite element method for 
parabolic problems. The paper goes far beyond what are explored in this paper, which deals 
exclusively with the h- vers ion of the DG method. In Section 3, a numerical example is given to 
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demonstrate the validity of our theory. Also included in Section 3 is the discussion of automatic 
time step control algorithm which is based upon the results presented in Section 2. The term 
||h 2 .D 2 u||/ n in (1.5) describes the spatial discretization error and it is of the second order due 
to the use of linear splines in defining the space V m as long as ||D 2 u(£)||/ n is well behaved. 
Unfortunately, for solutions of parabolic problems, this is not always the case. We investigate 
the behaviour of D 2 u and establish in Section 2 the spatial discretization scheme that achieves 
the quadratic convergence in space variable. 

2 Discontinuous Galerkin Method for Parabolic Problems 

Consider one dimensional parabolic problem (1.3). Even though the method will be outlined in 
one dimensional setting, extensions to higher dimensions are possible and straightforward. 

For 0 < a < 1 and q a nonnegative integer, define an index of singularity Q = For a 
positive integer N and T > 0, let 

£ = (^) Q > n = 0,1, . . . , N 

and 

tn = t' n T. (2.1) 

Define I n = (£ n _i, t n ], n = 1, 2, . . . , N, and let k n denote the length of I n so that 

fc n = [(%) Q -( r ^) Q }T, n = 1,2, ...N. 

Note that 

(2.2) 

by the mean value theorem. 

The solution u(x, t) of (1.1) is then approximated in t over each I n by a polynomial of degree 
q. For instance, with q = 1, let P\w denote the linear interpolatory projection of w € H§ = 
{n: D^v E L, 2 ,j = 0, 1, 2\v = OondQ.} in time onto Whk , viz, 

Piw(x,t) = — r — -w(x,t n - 1 ) + - — ^ n ~- w{x,t n ), for each t E / n . 

kn k n 

Note that P\ is bounded with respect to the norm || * ||oo,/ n where 

lk(f)||oo,/„ = maxlHOIIioo(O). 

tfc Jn 
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Since fl is assumed to be bounded, P\ is bounded with respect to |j * |j/ n also. Since Pi equals 
the identity on polynomials of degree < 1, examining the error terms in the Taylor expansion of 
order 0 and 1 respectively, we obtain, for each n = 1, 2, . . . , N, 

\\u- Piu\\ In < ( \\u t {t)\\ 2 dt, (2.3) 

J I n 

and 

\\u-Pxu\\ In <Ckl\\u tt \\i n . (2.4) 

For higher values of g, carrying out the Taylor expansion to the order g, the best possible 
estimate that can be achieved for the projection P q : Hq — > Whk is given by 

\W~ P q u\\ Jn \\u t {t)\\ 2 dt,Ckl +1 \\^^u\\ In }. (2.5) 

The type of estimate obtained in (2.5) plays a useful role in bounding the error in time 
variable for the DG finite element method for parabolic equations. The first term on the right 
of (2.5) is used to capture the error during the transient phase of the solution, while the second 
term is used for the solution later in time, upon which ||ut II 2 is expected to be smooth. The 
following lemma examines the first term on the right hand of (2.5). In [9], for integro-differential 
parabolic equations with memory, a similar singular behaviour in solution is treated. However, 
the analysis used in [9] restricts the construction of graded time partitions to constant and linear 
cases only. 


Lemma 2.1 Let 0 < a < 1, q a nonnegative integer and T > 0, we assume that t n , n = 1, . . . , iV 
are defined by (2.1). Then 

f S- a ds<C-^rr, 

J h ~ N « +1 

where C is a constant independent on N and for n > 1 ? 


max 

Kn<N 



*ds < C — , 
~ N 


where C is a constant independent on N. 


Proof: For n~ 1, 

i s ~“ ds = = o (( i)"«). 
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For 1 < n < N, we obtain 

h n s ~ a ds < fj n [{ I jtr) Q T]~ a ds as s~ a is decreasing over 7 n , 

= r -a (I ^l)-a<3[( n )Q _ ( n^l)Q] 

— — 2 ^') by the mean value theorem 

— c( -L)Q~ a Q sSll— — r i r (_n^\aQ n Q-i-aQ 

(n-i)“« L 'Jpr?T t n _!; n 

r* n q ^ N q r* 1 

“ ^ JV9+1 — ^ 7V9+1 — L JV • 

Q.E.D. 

Using Lemma 2.1 with n = 1 and 9=1, the first term on the right of inequality (2.5) is 
of order O(^), provided that ||u t (i)|| 2 = 0(t~ a ) for some a G (0, 1) and for t € I\. As stated 
earlier, it is reasonable to expect that ||u tt ||j n < oc for n > 1, so that (2.5) yields 

\\u-Piu\\ In =0(±), for each n = 1,... N. (2.6) 

Lemma 2.2 Let t n and k n be defined by (2.1). Then, with any positive integer N, 

(1 + log ^y-) 1 / 2 < y/2, for each n = 0, 1, . . . , TV. 

Proof: 

= (! +log 1_(^1) Q ) 1/2 

= (l-l0g(l-(^)«)) 1 /2 

<(l-l0g(l-(^i)<?)V2 

< (1 + (^)°) 1/2 < V2. 

The second to the last inequality is obtained from log(l — x) < — x for x < L 
Q.E.D. 

Lemma 2.2 will be used in Theorem 1.1 to guarantee the stability of the DG method currently 
proposed. 

For spatial discretization, let M be a positive integer and let x* for k — 0, 1, ... , 2 M denote 
the spatial knots which will be defined precisely later in this section. Denote by tpk(x) the linear 
spline over = [xjfe-i,Xfc+i] for k = 1, 2, . . . , 2 M — 1 and ipo(x) and <p 2 m(x) are linear splines 


7 



on [xo,xi] and [x 2 a/_ 1 ,x 2 m] respectively. Extensions to higher order splines are straightforward. 
At each time level t n , we approximate u(t n ,x) by 

2 M 

u n = u n (x) = U(t n , x) = £ £Vi(x), n = 0, 1, . . . , N, 

7-0 

where £? = ii 0 (xi). Now, equations ( 1 . 2 ) can be restated as follows: For n = 1, 2, . . . , AT, given 
U n ~ l find U = U\ In e P q {I n ) such that 

[ [(U t ,v) + a(U,v)}dt + (U n - 1 ’ + ,v n - 1 ’+)= f (f,v)dt + (U n - 1 -,v n ~ 1 ' + ) (2.7) 

J In J In 

for all v € P q (I n ) where C/ 0,_ = no- 

For demonstration, let us take a detour to see how equations in (2.7) take specific forms for 
q = 0 and q = 1 . First, consider the case g = 0 , ~i.e., constant in time. As £/ n = U n ~ = t/ n_1, + 
in this case, with U n = ]T? = 0 £p<pi(x), (2-7) reduces to 


£?[(¥><> ¥>j) - ¥>j)] = J2 C Vi) + (/> V’j). J = 0, 1, . . - , 2 


M 


( 2 . 8 ) 


7=0 


For q = 1, we let J7|/ n = <p n (x) + 


i=0 
t — t n — 1 


^ n (x) where cp„ = £;=o and = 


£?=o£i ViO*). As [/ n “b+ = ^ and U n ~ l '~ = <p n -i+$n-i, (2.7) becomes 


E?=o Vi) + kna&i, ¥>;)} + Zti tt' n {{VuVj) + W)} 

= L n (f(t), vj) + £i=o[ef’ n_1 + # ¥*), J = 0, 1, ■ • • , 2 M 

Ef=o Vi) + £i=i £•"{$(<*, ?;) + W )} 

= £ Jr (t - tn-l)(/(«). 3 = 0, 1, - • - , 2 M 


(2.9) 


The following theorem is a generalization of Theorem 1 . 1 . A proof can be obtained by 
applying minor changes in the proof of theorem 2.3, [ 2 ]. 

Theorem 2.3 Suppose that there is a constant 7 such that time steps k n satisfy k n < 7 fc n +i> 
n = 1 , . . . , N — 1 and let U n denote the solution of (1-2) approximating u at t n . Here u is 
approximated by a polynomial of degree q > 0 over each I n for n=l,...,iV — 1 . Then there is 
a constant C depending only on 7 and a constant (3, where px > (3hx and px is the diameter 
of the circle inscribed in K for all K G T, such that for n = 1,2,..., iV, 


u(t n ) - U n \\ 2 < C( 1 + log -^) 1/2 {max \\u - P q u\\j m + ||/i 2 D 2 u||/ n }. 


' k n 


m<n 


( 2 . 10 ) 
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We proceed to examine the assumptions of Theorem 2.3 relative to the current time and 
space discretization scheme. First of all. the condition k n < ~ik n +\ is satisfied with 7 = 1 as k n 
is increasing. Second, for one-dimensional problem (1.1), hx = px = \£Ik\ so that ^ = 1 for all 
fix- Hence the assumptions are fulfilled under the current time discretization scheme. Lemma 
2.2 guarantees the uniformly bounded property of (1 + log j^-) 1 / 2 , establishing the stability of 
the current DG method. 

By (2.5) along with Lemma 2.1, with any nonnegative integer q , we obtain, provided that 
\\u t {t )\\2 = 0(t~ a ) for some 0 < a < 1, 

||u-P 9 «|k=0(^). (2.11) 

It remains to examine the spatial discretization error term max*<* n \\h^ l D 2 u(t)\\j n in (2.10) for 
the current problem. If \\D 2 u(t)\\j n is bounded for each n and for each t. then using the graded 
time partitions described in Lemma 2.1, Theorem 2.3 transforms to the following: 


Theorem 2.4 Consider the parabolic problem (1.1) and assume that the initial value uq(x) is 
defined in such a way that ||u t (t)||2 = 0(£~ a ), for 0 < a < 1. Denote by U n the solution of 
(1.2) approximating u att n . Let time discretization {£„} be defined by (2.1) and (0, n) is divided 
into 2 m subintervals each of equal length. Also assume is bounded for each n. If 

q denotes the degree of polynomials used in approximating u in time variable, then for each 
n= 1,2, . . . , JV, 


In many practical problems, the assumption in the previous theorem that \\D 2 u(t)\\i n is 
bounded for each n and t may not be possible, -e.g., examine the solution of one dimensional 
problem in (1.3). This topic is now considered. 

For elliptic problem, 


—A u = f in $7 
u — 0 on T, 

it is well known (cf. see [6], p.92) that, for a smooth dQ, 


IMIi* 2 (n) ^ CII/II 2 * (2.12) 

Now for convenience, let us assume that / = 0 in (1.1). Then using the definitions of the 
respective norms and from (1.1) and (2.12), we obtain 

\\D 2 u(t)\\i n < llu(f)IMn) < C\\u t (t)\\ 2 . (2.13) 
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If ||i£f (^) [| 2 is bounded as t — * 0 + , then spatial discretization can be made arbitrarily with step 
size h and attain 

max/i 2 ||n t (f )||2 = 0(h 2 ), where h — 2~ M at each time level t n . 

t<t n 

In the case, ||iq (£)||2 — 0(t~ a ), (2.13) imphes that we must select a set of knots {£*.(£)} that 
depends upon t . The approach here is similar to the adaptive strategy employed by Eriksson 
and Johnson [2] in which the spatial increment h n is selected according to the size of ||u*(£)|| 00 ,/ n - 
Let h(t ) — max 1 < fc < 2 A i{xk(t) ~ x k-i(t)) for each * £ (0,T]. Then h(t ) is determined from the 
condition that h(t) 2 t~ a = 0(t 2 ) as t — > 0+. In terms of N, we require that h(t) 2 t~ a = 

■ q + 1 

For t = tj = (-fa) T, j = 1, . . . , AT, we have, for some C > 0, 

* 2 (<j)[(^)^n _a =cT (2.i4) 

Solving (2.14) for h(tj) with C = 1 for convenience, 

i a Q 

Finally, select M so that 

2 " > m <2 ' 16> 

for each time level tj , j — 1 and Xk(tj) = (note that M depends upon j) , for 

k = 1, . . . , 2 m — 1 to form the spatial partition points. Reflecting on what was just discussed, 
we finally have the following theorem which characterizes the convergence of the DG method 
both in time and in space variables in terms of N. 

Theorem 2.5 Consider the parabolic problem (1.1) and assume that the initial value uq(x) is 
defined in such a way that || ti-t (t) || 2 = 0(t~ a ), for 0 < a < 1. Denote by U n the solution of (1.2) 
approximating u at t n . Let time discretization {t n } be defined by (2.1) and (0, 7r) is divided into 
2 m subintervals each of equal length , where M is defined by (2.16). If q denotes the degree of 
polynomials used in approximating u in time variable , then for each n = 1, 2, . . . , N 1 

ll«(U-%|| 2 = 0( ] T T + ] ^). (2.17) 

Note that equality (2.15) implies that less partition points in space would be required as j — > N. 
This characterizes the situation described previously typical of parabolic problems in which 


10 


solutions become smoother after certain transient periods. Also note that, at each time level 
tj particularly for t near 0, it is possible to reduce computational cost by relaxing the size of 
spatial elements over the regions where u(^tj) is well behaved. This will be made more evident 
when a numerical example is discussed in the next section. 

3 A Numerical Example and Automatic Time Step Control 

In this section, the graded partitions described in the previous section is tested in the following 
standard one dimensional problem. Find u such that 

U t (x , t) - U XX (X , t) = 0, 0 < X < 7T, t £ i? + , 

u(0, £) = u(7r,t) = 0, t € i2 + , (3.1) 

u(x , 0) = 7 r — x, 0 < x < 7r 

The initial function is incompatible with the zero boundary conditions in (3.1). The exact 

solution is given in (1.3) and ||wt(t)||L 2 (n) = 0(i~5). Thus, the indices of singularity Q — 
are 4 and 8 respectively for constant and linear approximations in time variable. Schotzau and 
Schwab (example 7.1 [12]) discuss a similar situation in which the initial function is taken to be 
u(x,0) = 1. For this case, the grading function h is selected as h(t) = t 3 ( 2 ?+ 3 ) in [12]. With N 
denoting, as before, the number of partitions in time variable, Schotzau and Schwab’s approach 
results in t\ = (7^) 9 and ti = (*) 15 for the constant and linear approximations respectively in 
[12]. In the present method, as ||«t(^)llL 2 (n) = 0(t~%) with u(x,0) = 1, t\ — (^) 4 and t\ ~ (jf) 8 
for the constant and linear approximations. This shows that the present method provides more 
relaxed time increments at the beginning, which is critically important in order for numerical 
computation to proceed with a reasonable cost. 

We select T — 5 and N = 10 so that the solution is sampled at each time level tj = (^)^T, 
j = 1, . . . , 10. First, the case for constant approximation in time, q ~ 0, is considered. The exact 
solution and the approximate solutions are plotted in Appendix for t n , t n = 1, . . . , 9. For q — 0, 
the errors are dominated from the time discretization term which is of order O(jj). In Fig. 1.1 
with t\, from (2.15), h(t\) 2^ .0058. Hence we take M ~ 9 in (2.16) to guarantee the convergence 
rates in Theorem 2.5. It is interesting to note that, for the example under consideration, we can 
relax M significantly while maintaining an overall quality of the numerical solution. Also recall 
that the size of spatial discretization can be relaxed over the region away from x = 0 in this case 
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where u(x,ti) is well-behaved, but, of course, this requires a priori knowledge of u(x.t i). For 
q — 1, h(t\) ~ .0002 from (2.15) so that choose M = 14. In this case, (2.17) gives 

Nil) - Uih = O(dj). 

Table 1 below lists L 2 errors at each time level t ny n= 1, 2, . . . , 10. 



Table 1 

n 

q = 0 

9 = 1 

1 

.0133 

.0037 

2 

.0571 

.0020 

3 

.0673 

.0089 

4 

.0729 

.0087 

5 

.0763 

.0086 

6 

.0780 

.0081 

7 

.0878 

.0076 

8 

.0933 

.0072 

9 

.0702 

.0059 

10 

.0850 

.0073 


We close this paper by considering an automatic time step control scheme for 
for parabolic problems having non-smooth data. Johnson considered a similar 
step control method for a class of stiff ordinary differential equations in [7]. 
extended to an automatic time control method for a nonlinear parabolic problem in [1]. Let us 
recall the following algorithm in [1] for choosing the time step size k n : 


the DG method 
automatic time 
The result was 


Algorithm: [1] Given a tolerance 6 > 0, 

1. Choose k n = fc n -i? 


2. Given fc n , compute the corresponding approximate solution U(t n ). 

3. If 


;L< l)l|2<|, 

where 7 is a suitable constant, 7 ~ 2 or 3, then stop and accept the time step k n . Otherwise, 
decrease or increase k n by a factor of, say, 2 and return to step 2. 


12 



To demonstrate that the algorithm above may become expensive in the presence of a singular 
transient phase, consider the following: Suppose that. S = 10^ 5 and ||w f (<)|[ 2 = 0(t~ 1 ^ 2 ). Ap- 
proximating \\U(t\ ) - U (to) H 2 by and decreasing the time step size by a factor of 2 beginning 
with ki = it would take 8 steps to find k\ which produces the numerical solution within the 
prescribed tolerance. The current approach is to estimate the parameter a in ||xxt(t) || 2 = 0(t~ a ). 
The following algorithm contains a process of adjusting each time step (step 4), if necessary, as 
computation progresses, but, in principle, if a is estimated accurately in step 1, the amount of 
adjustments would be minimal. 


Algorithm: Given a tolerance 5 > 0, 

1. Choose t\ = k\ small. Compute || L7 (ti) — f/(0)||2. As 


so a is estimated by 


m*i)ii 2 * =< ckr , 


, Iniiu(ti)-D'(0)ii 2 
0 1 ' 


2. Time discretization error is controlled by requiring 

f \\u t (t)\\ 2 dt<6. 

Jh 

t\ is adjusted by 

I-' < «, 

1 — a* 

and subsequently choose the number N of time partitions from 


where q is the degree of polynomials used in time approximation. 

3. With Q* = £*., define 

U = (^) q ' t , i = 2,3, . . . N, 
to complete the partitions 0 < t\ < < * * ■ tN = T. 

4. For i > 2, compute U(U) and check an accuracy requirement 


(3.2) 


W{U)-U(U-i)h<8. 

If this inequality is satisfied, then proceed to compute U(ti+\). If not, adjust tj+i by 
decreasing ki + 1 = — ti by the factor of, say, 2. 


13 


References 


[1] K. Eriksson and C. Johnson, Error estimates and automatic time step control for nonlinear 
parabolic problems , /, SIAM Jl. Num. Anal., Vol. 24, No. 1, (1987), pp.12-23. 

[2] K. Eriksson and C. Johnson, Adaptive finite element methods for parabolic problems I: a 
linear model problem , SIAM Jl. Num. Anal., Vol. 28, No. 1, (1991), pp. 43-77. 

[3] K. Eriksson and C. Johnson, Adaptive finite element methods for parabolic problems II: 
Optimal error estimates in L 00 L 2 and L^L^, SIAM Jl. Num. Anal., Vol. 32, (1995), pp. 
706-740. 

[4] K. Eriksson and C. Johnson, Adaptive finite element methods for parabolic problems IV: 
Nonlinear problems, SIAM JL Num. Anal., Vol. 32, (1995), pp. 1729-1749. 

[5] C.M. Elliott and S. Larsson, Error estimates with smooth and nonsmooth data for a finite 
element method for the Cahn- Hilliard equation, Math. Comp., Vol. 58, (1992), pp. 603-630. 

[6] C. Johnson, Numerical Solution of Partial Differential Equations by the Finite Element 
Methods, Cambridge University Press, (1994). 

[7] C. Johnson, Error estimates and adaptive time-step control for a class of one-step methods 
for stiff ordinary differential equations, SIAM Jl. Num. Anal., Vol. 25, (1998), pp. 908-926. 

[8] H. Kaneko and Y. Xu, Superconvergence of the iterated Galerkin methods for Hammerstein 
equations, SIAM Jl. Num. Anal., Vol. 33, No. 3, (1996), pp. 1048-1064. 

[9] S. Larsson, V. Thomee and L.B. Wahlbin, Numerical solution of parabolic integro - 
differential equations by the discontinuous Galerkin method, Math. Comp., Vol. 67, (1998), 
pp. 45-71. 

[10] X. Makridakis and I. Babuska, On the stability of the discontinuous Galerkin method for 
the heat equation, SIAM Jl. Num. Anal., Vol. 34, (1997), pp. 389-401. 

[11] L. Schumaker, Spline Functions: Basic Theory, John Wiley & Sons, NY (1981). 

[12] D. Schotzau and C. Schwab, Time discretization of parabolic problems by the hp-version 
of the discontinuous galerkin finite dement methods, SIAM J. Numerical Anal., Vol. 38, 
(2000), 837-875. 


14 


[13] D. Schotzau and C. Schwab, An hp a priori error analysis of the DG time stepping methods 
for initial value problems, Calcolo, Vol. 37, (2000), pp. 207-232. 


15 



Solidline=exact: Dotline=approx Soiidline=exact: DoHine=approx Soiidline^exact: Dotline=approx 


Appendix 

Fig 1.1 Discontinuous Gaterkin Finite Element Method for Parabolic Problem - constant time 
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Fig 1.4 Discontinuous Galerkin Finite Element Method for Parabolic Problem - constant ti 


n=4 

Error* 0.0729 



Fig 1 .5 Discontinuous Galerkin Finite Element Method for Parabolic Problem — constant ti r 
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Fig 1 .6 Discontinuous Galerkin Finite Element Method for Parabolic Problem — constant tir 
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Fig 1 .7 Discontinuous Galerkin Finite Element Method for Parabolic Problem - constant time 
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Fig 1 .8 Discontinuous Galerkin Finite Element Method for Parabolic Problem — constant time 
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Fig 1 .9 Discontinuous Galerkin Finite Element Method for Parabolic Problem — constant time 
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